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We report the results of muon-spin relaxation measurements on the low-dimensional antiferro- 
magnet Rb4Cu(Mo04)3. No long-range magnetic order is observed down to 50 mK implying a ratio 
Tn/J < 0.005 (where J is the principal exchange strength along the spin chains) and an effective 
ratio of interchain to intrachain exchange of | J±/J| < 2 x 10"'', making the material an excellent 
realization of a one-dimensional quantum Heisenberg antiferromagnet. We probe the persistent spin 
excitations at low temperatures and find that ballistic spin transport dominates the excitations 
detected below 0.3 K. 

PACS numbers: 75.10.Pq, 76.75.+i 



The one-dimensional S — 1/2 quantum Heisenberg an- 
tiferromagnet (IDQHAF) described by the Hamiltonian 
H ~ JJ2 ^i ' '^i+i continues to be of much interest, not 
only as a system in which quantum fluctuations control 
the low-temperature physical behavior—, but also due to 
the question of whether the spin transport in this system 
is ballistic or diffusive^i^. However, the experimental ver- 
ification of theoretical work relies on the existence of ma- 
terial systems whose interactions are well approximated 
by the assumptions of the model. Such systems include 
the inorganic compound^i^ Sr2Cu03, the organometallic 
coordination polymer^ Cu(pyz)(N03)2 and the organic 
radical-ion sah^ DEOCC-TCNQF4. 

Recently^ the inorganic system Rb4Cu(Mo04)3 was 
identified as another candidate quasi-one-dimensional 
5 = 1/2 magnetic system, where the principle magnetic 
exchange occurs along one-dimensional (ID) chains, with 
exchange strength J. Real systems do not comprise com- 
pletely isolated chains and there are also interchain in- 
teractions described by an effective exchange coupling 
J±. Such a coupling, although weak in comparison to 
J in many materials, controls the thermodynamics and 
leads to long-range magnetic order (LRO) at a nonzero 
temperature T^. The ratio T^/ J is therefore a good in- 
dicator of the degree of isolation of a low-dimensional 
magnetic system. In fact, quantum Monte Carlo sim- 
ulations of this effect provide a means of estimating the 
effective interchain coupling J± in a quasi-lD magnet via 
the expression^ 



\Ji 



4c\/l^(M)+5lnln(M) 



(1) 



where c = 0.233 and A = 2.6. It is important to note that 
the ID nature of the spin system promotes strong ther- 
mal and quantum fluctuations which lower the magnitude 
of Tn, resulting in quite long-ranged correlations building 
up in the chains just above the ordering temperature, and 
also reduce the magnitude of the magnetic moment^. 



FIG. 1: The principal exchange paths in Rb4Cu(Mo04)3 
(viewed along the &-axis) showing strong exchange (with 
strength J) and weak exchange (strength J') between Cu^"*" 
ions through Cu-O-Mo-O-Cu paths. Weak exchange bonds 
between Cu^^ ions and apical oxygens are dotted. 



These factors often make the identification of magnetic 
ordering difficult using thermodynamic probes. We have 
shown previously^'!'^'!! that muons are often sensitive to 
LRO that can be very difficult to detect with conven- 
tional probes. At temperatures above Tn muon-spin re- 
laxation (/^+SR) may also be used to probe the nature 
of the spin excitations and their transport. In this paper 
we investigate the magnetism in Rb4Cu(Mo04)3 using 
implanted muons. We search for long range magnetic or- 
der down to 50 mK and investigate the spin excitations 
of the system. We argue that the system remains dis- 
ordered at all measured temperatures and the persistent 
mode that dominates the dynamics we measure below 
0.3 K is ballistic in character. 

The structure'^ of Rb4Cu(Mo04)3 is layered, with tri- 
angular lattice sheets oi S = 1/2 Cu^"'' ions in the a-c 
plane [Fig. [T], separated by layers of Rb+ and M0O4. 
Although this constrains the principal exchange to occur 
along paths in the a-c plane, the two dominant exchange 
paths, with strengths J and J', arc quite different. The 
strength of the exchange follows from the disposition of 
the Cu^+ dx2-y2 orbitalsS'. The local environment of the 
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FIG. 2: Example ZF ^+SR spectra measured at T = 0.05 K 
and r = 7 K. 



copper is a CuOs distorted square pyramid which causes 
the d^2_y2 orbitals to he in the base of these pyramids. 
Strong a bonds therefore hnk the dx2-y2 orbitals on the 
Cu^"*" and the Px,y orbitals on neighbouring 0^~ ions that 
form the basal plane of the pyramids, with one strong O- 
Cu-0 bond in the a-c plane and one directed along b. The 
main superexchange pathway, with strength J, therefore 
goes via two a bonds in the a-c plane, along identical zig- 
zag Cu-0-Mo-O-Cu paths, as shown in Fig.[TJ The other 
possible exchange path in this plane, which has strength 
J', is suppressed since it involves only one a bond and 
also the 0-Cu bond [shown with dotted lines in Fig. [T] 
that involves the apical oxygen of a CuOs unit, whose pz 
orbitals enjoy very little overlap with the Cu^^ dx2^y2 or- 
bital. In this way, the spin system appears to be formed 
from one-dimensional chains with coupling J, with weak 
interchain coupling J' in the a-c plane and coupling be- 
tween planes that is still weaker. This has led to the 
material being advanced as a strong candidate to show 
quasi-one-dimensional magnetic behavior. Magnetic sus- 
ceptibility and heat capacity measurements^^ confirm that 
the dominant magnetic exchange occurs along strongly 
interacting antiferromagnetic chains with an exchange 
constant J = 10.0 K. No signature of LRO was detected 
down to 0.1 K, with only broad humps in susceptibility 
and heat capacity in the 4 K region, typical of the build 
up of short range correlations above T^. 

In a muon-spin relaxation (/i+SR) experimenlji^, spin- 
polarized positive muons are stopped in a target sample, 
where the muon usually occupies an interstitial position 
in the crystal. The observed property in the experiment 
is the time evolution of the muon spin polarization, the 
behavior of which depends on the local magnetic field 
B at the muon site, and which is proportional to the 
positron asymmetry function A{t). Our /i+SR measure- 
ments were made using the MuSR spectrometer at the 
ISIS facility, UK. For the measurement a polycrystalline 
sample of Rb4Cu(Mo04)3 was mounted on a silver sam- 



ple plate using a small amount of vacuum grease and the 
plate was attached to the cold finger of a dilution refrig- 
erator. 

Example zero-field (ZF) /i+SR spectra are shown in 
Fig. [2] Across the entire measured temperature range 
(0.05 < T < 50 K) A{t) is found to relax exponentially 
with small relaxation rates of A < 0.03 MHz. Data were 
fitted to 



A{t) = Aoe-^' + Ab 



(2) 



where Abg is a small, constant offset reflecting the con- 
tribution of those muons that stop in the sample holder 
or cryostat tail. Oscillations, which would be charac- 
teristic of a quasi-static local magnetic field at the muon 
stopping site, are not observed at any temperature. How- 
ever, we do observe a cross-over in magnetic behaviour 
in the system upon cooling below Tco ~ 5 K. Below this 
temperature the initial asymmetry drops by an amount 
AA = 0.7(1)% (translating to a loss of roughly 3% of the 
total muon polarization) as shown in Fig. ^h) . Also evi- 
dent is an increase in the relaxation rate A, which begins 
to rise below 4 K, increases steadily and then flattens off 
below 0.3 K [Fig. i^a)]. 

These results are indicative of a transition or cross-over 
below Tco to a regime characterized by larger magnetic 
fields, a broader distribution of fields or (perhaps most 
likely) an increased correlation time. Assuming the fast- 
fiuctuation limit (justified below) we expectJ^ the relax- 
ation to vary as A ex A^r, where A is proportional to the 
second moment of the field distribution via A^ = lfi{BV), 
7^ is the muon gyromagnctic ratio, B^ is the magnetic 
field at the muon site and r is the correlation time. The 
small fraction of asymmetry missing below T^o therefore 
refiects a fraction of muons which are quickly depolarized 
within the first ^0.1 /xs and the increased relaxation of 
the remainder of the muons also indicates an increase in 
the magnitude or width of their field distribution via A 
or correlation time via t. The small size of the muon re- 
laxation rate (leading to a decay time that is rather slow 
in comparison to the length of our time window) means 
that reliably extracting a value of the asymmetry base- 
line is difficult. We do not, however, see any evidence for 
an increase in the apparent baseline on cooling through 
Tco as would be expected in a transition to static mag- 
netism in a polycrystalline sample. (Such an increase 
occurs since 1/3 of muon spin components would be ex- 
pected to lie along the direction of a static local field and 
then not be relaxed.) This suggests that the magnetic 
state is characterized by significant dynamics in the local 
field at the muon site(s) at all temperatures. 

To understand the nature of the muon relaxation, mea- 
surements were made in applied magnetic field, directed 
parallel to the initial muon spin polarization. The appli- 
cation of fields of i? < 3 mT leads to a repolarization of 
the asymmetry. In view of the form of the ZF spectra, 
fits were carried out to a dynamicized Kubo Toyabe (KT) 
function^"' with a field fiuctuation rate v = 1/r. Fits to 
this function in applied field with a second moment of 
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FIG. 3: (a) The relaxation rate A as a function of temper- 
ature is seen to increase upon cooling through Tco « 5 K. 
(b) The initial asymmetry A{t = 0) also shows a reduction 
upon cooling through Tco (dotted line), (c) The evolution of 
the relaxation rate A with applied field (above B — 3 mT) at 
T = 0.1 K. All fits are described in the main text. 



the field distribution A « 0.1 MHz and a fluctuation rate 
1/ Ki 0.7 MHz describe the data at all measured fields be- 
low 5 mT. The ratio v/A w 7 places the relaxation in the 
fast fiuctuation limit, justifying the assumption above. 

An obvious explanation for the change of behaviour 
around 5 K would be a transition to a state of LRO at Tco. 
Although one might expect that LRO would lead to os- 
cillations in the asymmetry, there are several cases where 
such oscillations are not observed at a magnetic transi- 
tion. Similarly the absence of an apparent increase in 
baseline does not preclude the occurrence of LRO either. 
However, we believe LRO unlikely. If Tco = 4.8(3) K 
[which we estimate from the behavior of A{t = 0)] is iden- 
tified with the magnetic ordering temperature Tn then we 
have a ratio T^/J = 0.48(3). The model in Eq. ^ refers 



to a chain coupled weakly in the two perpendicular in- 
terchain directions, assumed identical, with strength J±. 
Although the interchain coupling in Rb4Cu(Mo04)3 is 
more complex, by applying Eq. ([T]) we obtain a predic- 
tion for an effective interchain coupling |Jj.| = 3.4(2) K 
which should be expected to be of the same order of 
magnitude as J' identified above. Using this estimate we 
find \J±/J\ = 0.34(2), which would imply a rather three- 
dimensional magnetic material with little of the quite sig- 
nificant exchange anisotropy suggested by the structure 
or the previous magnetic measurements®. In addition, 
chain mean field theory gives u s a predicted staggered 
momenti^ mg = 1.017a/| J_l|/ J w 0.6/ub, which is likely 
to be detectable. It therefore seems unlikely that the 
weak effects we observe with muons are LRO. Another 
possibility is that small phases of impurities are present 
and that these are ordering. In this picture the small 
loss of asymmetry would be due to the muons which stop 
in these regions and that the increased relaxation effect 
would reflect the magnetic field of the impurity regions in 
the rest of the sample. This seems highly unlikely in that 
if such impurities are detectable with /^"''SR their volume 
fraction must exceed several percent and their presence 
would likely have been detected with x-ray diffraction. 

It is most probable that the effects we observe at 
Tco are due to the rapid growth of correlations in 
Rb4Cu(Mo04)3 around 4 K. We note that broad peaks in 
magnetic susceptibility and in specific heat are both ob- 
served in the 4 K region^. These are due to the build up 
of magnetic correlations within the chains and reflects the 
fact that the energetic cost of flipping a spin is of order J. 
This build up could conceivably lead to local regions of 
well correlated spins which give rise to an increased muon 
relaxation. If that is the case then we must conclude 
that no LRO is observed down to 50 mK, placing strong 
constraints on the effective interchain coupling and the 
size of the staggered moment. We find Tn/ J < 0.005, 
|Jj.| < 0.02, \J±/J\ < 2 X 10-3 and toq < 0.05AiB. For 
comparison, we note that Cu(pyz)(N03)2 has \J±/J\ = 
4.4 X 10-3, SraCuOa has \J±/J\ = 0.93 x lO'^ while 
DEOCC-TCNQF4 has \J±/J\ < 0.06 X 10-3. These ob- 
servations imply that Rb4Cu(Mo04)3 is a better example 
of a IDQHAF than Cu(pyz)(N03)2 and perhaps compa- 
rable to Sr2Cu03. 

That the physics of Rb4Cu(Mo04)3 is well described 
by a ID model of spin exchange with small interchain 
interactions allows us to relate our results to the spin ex- 
citations expected for such a model. Below Tco [Fig[31[a)] 
A increases as the temperature is reduced, and saturates 
at a constant value Asat = 0.025 MHz below 0.3 K. This 
trend of an increased relaxation rate and low-T persis- 
tent dynamics is often seen in the /i^SR of complex sys- 
tems with low-temperature dynamic a^^i^^ but its origin 
is not fully understood. The general behavior can arise 
because of the existence of two relaxation channels, one 
of which is strongly T dependent with a correlation time 
Ts(T) while the other shows little variation with T and 
has a correlation time r^ . In the presence of two compet- 



ing relaxation processes the one with the shorter corre- 
lation time wins out, giving the smaller relaxation rate. 
At higher temperature, therefore, we have Ts{T) ^ t„ 
which results in a strongly T-dependent relaxation. We 
may gain an insight into the temperature dependent part 
of A by by appealing to the results of previous measure- 
ments on highly ideal IDQHAFs. The muon is a local 
probe, so its response is due to a sum over fluctuations 
at all wavevector q. However, NMRf^ of Sr2Cu03 and 
//+SRl of DEOCC-TCNQF4 revealed that the low tem- 
perature relaxation was dominated by spinon excitations 
close to g = tt/o where low T corrections to the theory 
of spin correlations for the S — 1/2 IDQHAF result in a 
rise in A that follows A ^ y/\n{4~5J/T) . Such a tempera- 
ture dependence provides a good description of the rapid 
increase in A at low temperatures in Rb4Cu(Mo04)3 as 
shown in Fig. El^a)^^. The increase in relaxation rate 
arises from the lengthening of the correlation time Ts 
with decreasing temperature. This increase continues un- 
til Ts becomes much longer than the correlation time r„ 
of the T-independent relaxation channel with persistent 
dynamics, which then dominates the relaxation at the 
lowest temperatures, resulting in a relaxation rate Asat- 

Finally we turn to the nature of the propagation of the 
persistent spin fluctuation in Rb4Cu(Mo04)3. Whether 
the excitations in IDQHAFs are ballistic or diffusive re- 
mains a key question in the physics of spin chains and one 
that can be addressed with our LF ^+SR measurements^. 
This is because the spectral density for ballistic motion 
is expected to vary with frequency uj as fij-j) ^ \n(J/uj) 
while ID diffusive transport results in a function /(w) ~ 
uj-i. In applied fields B > 3 mT at T = 0.1 K we 
begin to repolarize the missing fraction of asymmetry 



and obtain a relaxation that persists up to the highest 
applied field of 250 mT. The result of fitting this per- 
sistent relaxation rate to an exponential relaxation func- 
tion as a function of applied field is shown in Fig. ^c) . 
We see that the diffusive model, with A oc B^^-^ pro- 
vides a poor description of the data. In fact, fitting 
to an unconstrained power law form A oc i?" results in 
n = 0.20(2), well below the theoretical value of n = 0.5 
and the slightly reduced value of n = 0.35 found^ in 
DEOCC-TCNQF4. Turning now to the model of balhstic 
transport, we note that the muon can couple to electronic 
spin density via both dipolar and hyperfine interactions, 
and for ballistic spin transport we expect^^ a relaxation 
rate X{B) = ^ [3^VK) + (5^' + 7^')/(^e)] , where 
■^/j = 7/j^7 '-^c = "fcB, 7o is the electron gyromagnetic 
ratio and A and D are the scalar and dipolar hyperfine 
coupling parameters. As shown in Fig. |3l[c), a good fit 
to the data is achieved by fitting A oc ln(J/wo), which 
assumes that the second term in the expression for X{B) 
dominates. Our inference, therefore, is that the modes 
which cause the saturated relaxation at low temperature 
show behavior suggestive of ballistic spin transport. 

In conclusion, our results demonstrate that 
Rb4Cu(Mo04)3 is a good example of a quasi-lDQHAF 
system, with a low ratio of exchange strengths | J_l/J|, 
putting it amongst the best realizations of this model. 
The spin excitations are suggestive of a contribution 
from modes around q — tt/o in the region 0.3 < T < 4 K 
and a persistent, low temperature contribution which 
has a ballistic character. 

Part of this work was carried out at the STFC ISIS 
facility and we are grateful for the provision of beamtime. 
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